Computation of Resistive Wakefields 



o 
o 

(N 
G 



43 ■ 

Oh: 

I 

O 1 

o : 

: 

1/3 1 

O ' 
• I— I ■ 

C/3 ■ 

43 : 

Ph. 



> 

o 

OO 

q 

o 
o> 
o 



X 



Adina Toader and Roger BarlowS 
School of Physics and Astronomy, The University of Manchester, Manchester M13 9PL, United Kingdom 

We evaluate longitudinal resistive wakefields for cylindrical beam pipes numerically and compare 
the results with existing approximate formulae. We consider an ultra-relativistic bunch traversing a 
cylindrical, metallic tube for a model in which the wall conductivity is taken to be first independent 
and second dependent on frequency, and we show how these can be included simply and efficiently 
in particle tracking simulations. We also extend this to higher order modes, and to the transverse 
wakes. This full treatment can be necessary in the design of modern nano-beam accelerators. 

PACS numbers: 



I. INTRODUCTION 

The requirements of modern accelerators are placing 
increasing demands on technology, with small precisely- 
defined beam bunches passing through physically small 
beam pipes and collimators. This means that the effects 
of wakefields produced by induced charges and currents 
are becoming increasingly important and require accu- 
rate calculation to ensure they do not dilute the emit- 
tance of the beam. 

The wakefields from particles in a bunch may affect 
other particles in the same bunch (intrabunch wakes) and 
also particles in subsequent bunches (intcrbunch wakes). 
For physically small apertures the effects of short-range 
intra bunch forces become important and need to be con- 
sidered, as well as the more familiar long range inter- 
bunch effects. 

For a relativistic particle in a perfectly conducting uni- 
form beam pipe the wake is zero. For a non-uniform 
beam pipe of finite conductivity the effect is thus sepa- 
rated into geometric and resistive wakes, the second of 
which is considered here. 

The wake is due to the Lorentz force F — q(E + vAB), 
however for the longitudinal component fn the second 
term is zero and only En need be considered. For a purely 

resistive wake, the wake field, E+ vAB, is constant and 
the integrated effect of the complete transit of the bunch 
through a section of pipe, generally called the wake po- 
tential, is obtained just by multiplying the wake field by 
the pipe length. The wake field is found [l| by writing 
down Maxwell's Equations in the beam aperture and in 
the beam pipe and matching them subject to appropriate 
boundary conditions. These solutions can be written as 
a sum over angular modes. For many purposes only the 
knowledge of the leading modes (m = and m = 1) is 
adequate, however as requirements become more strin- 
gent a technique is needed that will work in the general 



Finding solutions of Maxwell's equations is more read- 
ily done in frequency (or wavenumber) space rather than 
physical space, as differentiation becomes multiplication. 
The Fourier Transform of the wake is the impedance 



F(k) = / F{s)e lks ds F{s) = 



1 

2^ 



F(k)e 



*Also at The Cockcroft Institute, Daresbury 
Science and Innovation Campus, Warrington, 
WA4 4AD, United Kingdom; Electronic address: 
Adina.Toader@manchester.ac.uk, Roger.Barlow@manchester.ac.uk 
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where s is the distance along the beam axis between a 
leading particle which creates the field and a trailing (wit- 
ness) particle which feels the effect. 

For many purposes a knowledge of the impedance suf- 
fices. The 'kick factor' averaged over all the bunch is also 
often useful and expressions are given in the literature. 
However the kick factor only gives the mean effect for the 
whole bunch and in particle tracking codes an expression 
is needed for the physical impulse that one leading parti- 
cle has on another trailing particle. Existing approaches 
use expressions in different regimes where the divi- 
sion is somewhat arbitrary. We have found an approach 
which unifies, simplifies and speeds up the calculations, 
and makes the underlying physics clear. It includes short 
range and long range wakes with no artificial division be- 
tween them. It can also be extended to arbitrary angular 
order. 

Chao [H gives a general formula for the impedance, 
from which he obtains an expression for the physical 
wake which is valid in the long-range limit. Bane and 
Sands 0,0] extend this to shorter range though still mak- 
ing approximations. 

Gluckstern, van Zeits and ZotterQ consider resistive 
wakefields for circular, elliptical and rectangular pipes, 
but they only evaluate the impedance, not the physical 
wake, and only consider the lowest (m=0 longitudinal, 
m = 1 transverse) modes. Yokoya @ also considers beam 
pipes of general cross section, though his results are self- 
admittedly complicated and not easy to implement. Lut- 
man, Vescovi and Craievich Q generalise the results to 
elliptical beam pipes, including circular and planar aper- 
tures as special cases. Although a very general approach, 
they consider only leading modes and the examples they 
give are specific and not directly applicable to simulation 
programs. 
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II. THE LONGITUDINAL WAKE FOR m = 

The Fourier Transform of the m = mode of the lon- 
gitudinal component of the wakefield is given by 



E z (k) 



2q 



1 



2 



-(i + x)(i + 2is) 



where 



A(fc) 
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i + sgn(k)) 



(2) 



(3) 



with q the charge of the particle, b the radius of the tube, 
a the conductivity of the pipe, and c the speed of light. 
This assumes axial symmetry, the validity of Ohm's law, 
relativistic particles, and that the skindcpth is smaller 
than both the thickness of the pipe and the tube radius, 
but is otherwise general. This enables the Bessel function 

solution of F(k) to be replaced by a sinusoidal form, i.e. 
the asymptotic form of Jq(x) cx i~ 1,,2 e", as suggested 
by Chao([l(,p.43), which leads to the (2*5) term in the 
Equation [21 

It is convenient to introduce so, the scaling length 



so 



cb 2 

27R7 

and thus the dimensionless quantities K, the scaled 
wave number and s' , the scaled length 

K = Sgk 

, = s_ 

so 

It is useful to note that 



(4) 



(5) 
(6) 
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where the upper sign applies for positive k, the lower for 
negative k. 

To find the corresponding wakefield requires the in- 
verse Fourier transform. This has been much studied in 
the literature [l], 0], using various approximations forms 
of Equation [2] and evaluating them by a contour integral. 
By contrast we will do the integration numerically, en- 
abling us to provide a general technique without making 
approximations as has been done previously P, 0, H[ and 
to evaluate the different regions wherein such approxima- 
tions are valid. 

The back transform can be written 



E z (s) 
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E z {s) 



so 
2tt 



\Re[f even {k)] cos(ks) 
Im{fodd(k)]sm(ksf)dk (8) 
Re[f even {K)]cos(Ks') 
+Im[f odd (K)]sm(Ks'))dK (9) 



where the functions are the even and odd parts of E z (k). 
This separation is necessary to avoid problems with the 
alternating signs and the modulus operations in Equa- 
tions El 



A. The first order approximation 

In the limit of large b compared to sq and neglecting 
high and low frequencies Equation[2]can be approximated 

by 



The Fourier Transform is given by Chao Q as 



1 c - 



(10) 



(11) 



However we can also evaluate it numerically, in prepara- 
tion for more complicated forms of E z (k). We have 
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feven(K) = -\f{K)+f{-K) 

fo dd {K) = \[f(K)- f(-K) 



= -VK, 

= iV~K (12) 



where the functions are the even and odd parts of E z (k), 
having taken out a common factor of The even part 
is purely real and the odd part is purely imaginary, as 
they must be. 

A function like \/Ksm(K) presents problems for nu- 
merical integrals as the function is oscillating with 
increasing amplitude, and any summation technique 
(trapezoidal rule, Simpson's rule or Gaussian quadra- 
ture) will not work. However we know that the integral 
exists because it can be done analytically. We can per- 
form the K integral by first integrating with respect to 
s' , thereby dividing by a factor of K so that the oscilla- 
tions decrease in size and the numerical integration can 
succeed 



rs/s 

E z {x)dx = «o / E z (x')dx / 
Jo 

M [°° ( Re [feyenjK )] „■ 
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sin [Ks') 



cos (Ks 1 ) ) dK (13) 



Im[f odd {K )] 
K 



The result can then be differentiated numerically, us- 
ing an intermediate interpolating function, to give the 
desired wakefield. 

This can be verified and results are shown in Fig. [Tl for 
6=1 cm and a = 5.8 x 10 7 (Sim) -1 . These results are 
indistinguishable from those obtained by plotting Chao's 
formula, Equation II II 
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FIG. 1: The long range wake given by Equation [13] and Equa- 
tion [11] Q. 



B. A more accurate formula 

If the large s requirement is relaxed then a second term 
has to be included in the denominator and Equation [5] is 
approximated by 



E z (k) 



2q 



1 



(14) 



The even and odd parts are 



feven (-^ ) 



fodd(K) = - 



(2 }k) + k 
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(15) 
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E z (s, b) = ^§fU(s/s ). The result is shown in Fig.H to- 
gether with the lowest order approximation. They differ 
greatly, although above s' ~ 6 the agreement is actually 
quite good and Equation [TT]ia a good approximation for 
long range wakes. This transform can also be performed 
analytically, using contour integration as was done by 
Bane and Sands [3J and our Fig. [2] agrees with their re- 
sults, including the key value U (0) = —4. 



The full formula 



For the full version of Equation [2] one gets 



feven (-^0 



fodd(K) 



-16i 



(16) 



where we have introduced the dimensionless quantity 
£ = Sg/6 2 . Although this is no longer a universal curve, 
it can still be expressed as a function of two variables 
(s' and £) rather than the full set of three. The earlier 
approximation corresponds to the function at £ = 0. 
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FIG. 3: The exact longitudinal range wake for various values 
off. 



FIG. 2: The longitudinal wake in the second approximation 
(blue line) with the long range form (red line) for comparison. 

Although the wake is a function of three parameters 
(s, b and a), the use of the scaling length s enables 
it to be written as a universal function U(s'), where 



Fig. [3] shows how the function when evaluated with 
our numerical technique changes for different values of £. 
It can be seen that for values of below about 0.1 the ap- 
proximation is very good. For a copper beam pipe with a 
radius of 1 cm the scaling length is of order 20 microns, so 
£ is very small in all practical cases at present. However, 
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for possible future collimators with very low conductiv- 
ity and small radius it might need to be considered. It is 
also important to check that the simpler formula is valid 
before using it. 



Results are shown in Fig. [4] It can be seen that there are 
quite large differences in behaviour, even though only a 
few terms in the formula contain an m. The value at 
s = is given by U m {0) = -4(m + 1). 



III. LONGITUDINAL: HIGHER ORDER 
MODES 

For higher modes, using the same technique of match- 
ing the solutions of Maxwell's equations in the aperture 
and the pipe, and using the same approximations as for 
Equation [141 for m > 



2T 



1 



^2m+l ikb _ A 
m+1 k 



(17) 



where I rn is the charge moment of order m: 
J J p(r, 9)r m cos(md) dr dO. Any angular distribution of 
charges at a particular radius can be described in terms 
of these moments. Equation [T7] can be separated into 
odd and even parts, taking out a factor of I m s /b 2m+2 
as before for simplicity 
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FIG. 4: The longitudinal wake for various modes. 



The full formula for higher modes, corresponding to 
Equation [2 is 



4J„ 



fr2m+l ikb ( 2k I A^\ f-\i i \ im 
m + 1 y \ k J y 2A6 J kb 

(Note that Equation [TS] is valid only for m > 0.) 



(19) 



This can be separated into odd and even parts as for Equation [TB] 

feven (-^0 — 
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fodd(K) = -8i- 



ev / if-^(e + 2v^) + 2(^ T -e 



2 / s 2 



A 



(20) 



We show the dependence on £ in Fig. [5] for m = 1 
and m — 5. It can be seen that the dependence on £ 
increases for higher modes. However even for the m = 5 
mode shown the values of £ that would correspond to any 
reasonable beam pipe diameter and conductivity are still 
so small that the deviation is probably unimportant. 

It is frequently stated that the m = 1 wakefield is re- 
lated to the m = wake by a factor p-. From Fig.[4]it can 
be seen that the first two angular modes do indeed have 
the same shape, although higher modes are significantly 
different. The equality E 1 = -jpE is true for Equations 
RBI and [TTJ as can be see by making the substitution. 

However this is not true for more general formulae, 
Equations [T9land [TBI as is shown in Fig.[BJ For £ = 0.5, 
which is admittedly large, there is a clear difference be- 



tween the two curves. 



IV. TRANSVERSE WAKES 

Transverse wake effects are generally of more impor- 
tance than longitudinal wakes, especially for collimator 
studies. The transverse wakefield experienced by a parti- 
cle with transverse position r due to another particle at 
r' can be written as a sum over angular modes 

F T {r,6,s) =J2 rm ~ lr ' m ( fcos ( me ) - 6sin{m0))W^(s) 

m 

(21) 
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FIG. 5: The m = 1 (top) and m = 5 (bottom) wake as a 
function of a' and £. 




m= 
m= 1 



where 9 is the angle between the two particles and s is 
the distance between them. The property [3T] applies to 
the force components and not to the electromagnetic field 
components([l|, p. 56), and for the transverse wakes the 
magnetic part has to be included. 
The Panofsky-Wenzel theorem 



dF T 

dz 



applies term by term giving 



(22) 



(23) 



Thus the transverse wake at any order can be obtained 
by integrating the longitudinal wake. This is especially 
convenient for our method as this integral is what is cal- 
culated in the inverse transform. 



U(s') 




FIG. 7: The Transverse wake in the first and second approx- 
imation (red and blue curves respectively). 




m= 1 
m = 2 
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FIG. 6: The m = longitudinal wake, doubled, and the 
m — 1 wake for £ = 0.5. 



FIG. 8: Transverse wakes - various modes with £=0. 



The m = 1 transverse wake is shown in Fig. [7] for the 
first and second approximations (equivalent to Chao's 
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FIG. 9: Transverse (top) m = 1 and (bottom) m = 5 wakes 
with £ variable. 



formula and Bane and Sands' respectively). With our 
formula we can study the general case. Fig. [8] shows the 
different transverse modes, for £=0, and Fig.[9]shows the 
dependence on £ for m = 1 and m = 5. The dominant 



U(s') 

2 - 



m = 
— m= 1 



— ro 




FIG. 10: The m = transverse wake, doubled, and the m = 1 
wake for £ = 0.5. 



longitudinal effect is the m = mode, and the dominant 
transverse mode is for m = 1. There is, as Chao makes 
clear ([l[, p81) no linkage between modes of different or- 
ders to be obtained from the Panofsky-Wenzel theorem. 
However the m = and m = 1 longitudinal modes are 
connected by a simple factor of p- as discussed in Sec- 
tion (TTTJ In the literature (e.g. Q Equation 3) one finds 
the expression 



(24) 



This is exact for resistive wakefields only in the long wave- 
length approximation. We investigate the validity of this 
widely-used approximation in the general case. This dif- 
ference is small, but it is rather larger for the correspond- 
ing transverse wakes, as shown in Fig. [TUl 



V. AC CONDUCTIVITY 

In the classical Drude model for the ac conductivity 
case, in units of normalised wavenumber K — sok, we 
have 



1-iKT 



(25) 



where we have introduced the dimensionless relaxation 
factor r = ct/sq. For a 1 cm radius copper tube at room 
temperature the relaxation time is r = 2.7 x 10~ 14 s or 
ct = 8.1 fim giving r = 0.4, so we explore T values in the 
range to 5. In order to calculate the ac resistive wall 
wakefield the dc conductivity a in Equation 0] is replaced 
by the general form of the ac conductivity a. Hence, 



K\(l+K 2 T 2 )- 1 ^[iVT+t±VT zr i} (26) 




t±Wl + t 



with 



t 



\K\T 



Vi + K 2 r 2 



(27) 



(28) 



Following the same procedure employed on Section |TT] 
we express the inverse Fourier Transform of the m=0 
mode of the longitudinal component of the wakefield in 
terms of the even and odd parts. The simplest of cases 
given by Equation [12] becomes, 



j even (-^0 



I ' 



fodd(K) 



i 



where 



/ = 



\K\(1 + IPT 2 ) 1 /* 



(29) 



(30) 
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FIG. 11: The Long range wake for various values of F. 



The results are shown in Fig. [TT] for various values of 
r. If we consider the more accurate formula given by 
Equation 1141 our technique gives 



feven (-^0 



21^/T^t 



fodd(K) 



2i(§-l^/T+t 



(f-vr+t) +i 2 {\-t) 



(31) 



For the full version of the dc regime Equation [2] the 
odd and even parts have counterparts in the ac regime 
as 



feven (-^0 



-8 l 2 t(£ 2 y/l-t 2 + 2ilK^/T~t 
+Al 3 Ky/l~t)/d{K) 



fodd(K) = -8i l 2 t^ 2 t - 2£l(l - KVT+t) 
-4l 3 KVT+t + 2l 2 K 2 



/d(K) (32) 



with the denominator given by 
d(K ) = Vl-t 2 + 2£lKVl~^t + Al 3 KVT^t 
+ \ft - 2£l(l - KVTTt) - 4l 3 KVlTt + 2l 2 K 2 



For £ = one can see that Equationl3"2"lreduces to Bane 
and Sands' approximation, Equation [3T] as seen before. 
Fig. [T2l shows the results for various values of V. These 
wakes show a fairly strong dependence of T and the effects 
of AC conductivity may be important in a particular case. 
However £ will always be very small so we include it only 
for completeness and in practice one can probably use 
£ = as a good approximation. 

For higher order modes the impedance is given by 
Equation 1191 In the ac regime the odd and even parts 




FIG. 12: The Longitudinal wake in the second approximation 
for various F values. 



f even {K) = -8 l 2 t (£ 2 y/l - t 2 + 2£lKVT~t 
+Al 3 KVT~t)/d{K) 

fodd(K) = -8t l 2 t\£_ 2 t- 2^1(1 - KVTTt:) 



-4l 3 KVTTt + 4:l 2 ( 



K 2 



771+1 



(33) 
/d(K) 



with the denominator given by 
d(K) = (fVl - t 2 + 2£lKVl~t + 4l 3 KVl~t 



+4Z 



Ct - 2£l(l - KVT+t) - Al 3 KVl + t 

K 2 n -is 



TO + 1 



m£)K 



In Fig. [T3] we show the different modes dependance for 
two particular values of T, 0.4 and 5.0 and £ = 0. 




The r dependence of two modes, mode to = 1 and 
rn = 5 is shown in Fig. [14] also for £ = 0. The equivalent 
transverse wake is shown in Fig. [T5l It is clear that for all 
modes, increasing T increases the size of the wake effect, 
both longitudinal and transverse over the same distance. 
For the higher modes the effect is stronger. 



VI. IMPLEMENTATION 

The integrals used to generate the plots in this pa- 
per were performed using Mathematica 0]. Tables with 
6 x 1001 x 21 elements were written to file, covering the 
ranges < V < 5, < s' < 100 and < £ < 2 re- 
spectively. A separate table was used for each mode, 
and tables were written separately for longitudinal and 
transverse wakes. 

For ease of use we have defined a small object 
collimatortable written in C++ but portable be- 
tween Merlin, Placet, and hopefully other simulation 
codes. When created by a call to the constructor 
collimatortable (filename , Gamma, xi) it will read 
the full 3 dimensional table from filename and construct 
a one dimensional table, using parabolic interpolation be- 



30 

FIG. 14: The m=l (top) and m=5 (bottom) longitudinal 
wake as a function of V. 



tween the 9 closest points, for s' at this value of T and 
of £, both of which default to zero. It does this because 
when tracking particles through a collimator the values of 
s (or s') are different for each pair of particles (or slices), 
whereas the values of T and £ are the same, as they de- 
pend on the properties of the collimator, not the bunch. 
Then the one dimensional table can be used by the mem- 
ber function collimatortable: : interpolate (sprime) 
to find the value for any value of s' in the range. 

These data and program files are obtainable from the 
authors [Icj . 

A. MERLIN 

The MERLIN program [ll[ contains a general wake 
formalism which has been extended [l2| for high or- 
der modes (though only for axially symmetric aper- 
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FIG. 15: The m=l (top) and m=5 (bottom) transversal wake 
for various values of F, all for £=0. 



tuxes) and it is easy to incorporate these wake- 
fields. A class ResistivePotentials which inher- 
its from SpoilerWakePotentials calculates T and 
£ from the dimensions and properties given and 
uses collimatortable to read the files; functions 
Wtrans(z,m) and Wlong(z,m) use the interpolate 
function to get the wake experienced by the particles. 

For efficiency, the bunch is divided longitudinally into 
a number of slices N s (typically 100) and each wake func- 
tion is called only N S (N S — l)/2 times. 



B. PLACET 

The PLACET program Q contains only the lowest 
order m = 1 transverse mode, though it does include 
the effect of rectangular apertures using as an Ansatz 
the Yokoya factors y e ff = 0.822z/i + 0.411?/2, where y\ 
and j/2 are the displacements of the leading and trailing 
particle. y e ff is then used in the formula for circular 
apertures. 

In the current implementation a complicated system 
switches between different formulae depending on the 
sizes of the pipe and bunch. The integral resulting from 
the analytic transform of Equation [TJ] is calculated nu- 
merically. We replace these with a single call to the 
interpolate function. 100 lines of code in the original 



are replaced by 21 lines in the new version. The wake 
function is called for each slice for each macroparticle, 
N s N p /2 times. 

VII. EXAMPLES 

A. Simple collimator 

As a numerical example we consider aim long colli- 
mator made of Titanium, with conductivity a = 2.33 x 
10 6 (Om) — 1 . The radius is 1.4 mm and the relaxation 
time is 2.7 x 10~ 14 s. This gives s = 0.165 ^m, T = 0.49 
and £ = 0.00014. 
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FIG. 16: Merlin prediction for the deflection as a function of 
offset. 
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FIG. 17: Merlin prediction for emittance as a function of 
offset. 



This collimator and the beam bunch properties are 
those of the recent tests at SLAC End Station A [l3T ]. 
with X = 128.2 m, f3 y = 11.9 m, e x = 5.49 x 10" 9 m, 
e y = 3.44 x 10~ 10 m, initial beam energy p = 28.5 GeV 
and bunch length a z = 0.3 mm, with the sole difference 
of e x = € y = 3.44 x 10 -10 m when using Merlin. The 
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bunch of N p — 10 10 particles was modelled by 50000 
macroparticles in both simulations. 

We show the predictions of Merlin (PLACET) for the 
mean kick of the bunch in Fig. [T|5] (Fig. [T5)) and for the 
increased geometric emittance in Fig. [17] (Fig. [T9|) , re- 
spectively. The total emittance increase leading to the 
loss of luminosity at the IP caused by wakefields is due 
in part to jitter amplification and in part to the increased 
size of the bunch itself [l4[ • This latter effect is what we 
show. If one considers a bunch as a series of slices it is 
clear that different slices receive different angular kicks, 
the first slice receiving no kick at all, and this difference 
between slices increases the angular size of the bunch as 
a whole. 
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FIG. 18: PLACET prediction for the deflection as a function 
of offset. 



increase of emittance is 0.4 times coherent term produced 
by jitter. Here we show that it can be calculated explic- 

itly 

The results from Merlin and PLACET agree surpris- 
ingly well, given that one uses the Yokoya terms and the 
other does not. The overall shift in the emittance of 
Figs. IT71 and \W\ is due to the different random bunches 
and using only 50000 macro particles. Merlin results 
show that higher modes become important (only) for off- 
set > 6/2. 



B. Extreme collimator 

We consider aim long collimator made of graphite, 
with conductivity a = 7.14 x 10 4 (rim) _1 and very nar- 
row pipe of radius 0.125 mm. This gives Sq — 0.105 /mi, 
r = 2.85 and £ = 0.007. For comparison purposes we 
consider a beam bunch of the same properties as for the 
previous collimator, but with the emittance ten times 
smaller, while the displacements of the beam from the 
central axis are between and 0.115 mm. 
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FIG. 19: PLACET prediction for emittance as a function of 
offset. 



The deflection of Figs. [TH] and [TH] compounded with 
transverse jitter, may also contribute to the loss of lumi- 
nosity by an amount that depends on the betatron phases 
at the collimator and at the IP. In previous analyses, such 
as 11511. Ildl it was assumed that the incoherent intrinsic 



FIG. 20: Merlin prediction for the deflection as a function of 
offset. 
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FIG. 21: 

offset. 



Merlin prediction for emittance as a function of 
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for accelerators such as CLIC and the LHC. 
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FIG. 22: PLACET prediction for the deflection as a function 
of offset. 



Extreme Collimator 



-0.10 -0.05 0.00 0.05 0.10 

Bunch Offset (mm) 



FIG. 23: PLACET prediction for emittance as a function of 
offset. 



Again, the agreement is good (Figs. [20l - [23)) . Merlin 
shows that higher modes become significant sooner than 
they do for the previous case. Also it can be seen that 
more modes are needed for convergence when the bunch 
displacement is large compared to pipe radius. 



VIII. CONCLUSIONS 



Longitudinal and transverse resistive wakes can be cal- 
culated in simulations programs with a technique appli- 
cable for any general particle separation using pretabu- 
lated numerical Fourier Transforms and implemented in 
a way which is computationally efficient. Higher order 
angular modes can easily be included, and so can be AC 
conductivity effects. We have shown that this can be used 
by the Merlin, PLACET programs and other codes can 
be added in due curse. Results are given for two types 
of collimator, and show that the effects are as expected. 
This can be used in simulations of full beamline lattice 
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